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Abstract 

A generating function for products of Hermite polynomials is used 
to significantly simplify the evaluation of the integrals G n (p,p') oc- 
curring in the Gaussian model of multiple particle production. These 
integrals are crucial for studies of multiparticle effects in Bose-Einstein 
correlations. 
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1. Bose-Einstein correlations in multiple particle production processes 
are a field of intense study. For recent reviews see An important 

subfield is the study of multiparticle effects on multiplicity and momenta 
distributions. Here exactly soluble models are a convenient starting point. 
A soluble model, which is particularly popular, is the Gaussian model - 
13]. In this model all the integrations are Gaussian and consequently can 



be performed analytically. The generic integral is 

G n (p,p') = J pip, k 1 )dk 1 p(k 2 ,k 3 )---dk n -ip{k n -i ) i/), 
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where 



with 



P(p,p') 



V2nA' 



: exp 



2A 2 



-R 2 q 2 



(2) 



K 



p + p' 



p — p 



(3) 



The constants R 2 and A 2 are real and positive. For simplicity we start with 
the one-dimensional case. In the Gaussian model the transition to three 
dimensions is trivial. It is seen that the integrals G n are Gaussian and conse- 
quently can be performed analytically. One needs, however, a closed formula 
valid for any n. This is more difficult to obtain. Zimanyi and Csorgo || [7J 
following some suggestions from Pratt || constructed a system of recurrence 
relations for three series of parameters and after a rather lengthy calculation 
got a formula for the integrals, which can be written in a fairly simple form 
]l2fl . Wiedemann [|13| applied successfully the MATHEMATICA program di- 
rectly to the integrals G n . In the present note we point out that a very short 
derivation of the formula for the integrals G n can be given, if a mathematical 
identity for Hermite polynomials is made use of. This work is an extension 
of the method used in [§, §]. 

2. An identity, valid for \z\ < |, is 



00 z n 1 
&(x, y, z) = —H n (x)H n (y) = ■. exp 

(4) 

A convenient starting point to check this formula is the definition of the 
Hermite polynomials in the form p4| f] 

2 n r+°° ,2 
H n (z) = — (z + it) n e- 1 dt. (5) 

y7T J-oo 

After substituting it into ([|) the summation over n gives an exponential and 
it is enough to perform a double Gaussian integral to get the right-hand side. 
The function y, z) is defined by this formula. 

3. For any orthonormal set of functions ifj n {p) and any set of nonnegative 
numbers A n satisfying the normalization condition 



-Az 



z(x 2 + y 2 ) — xy 
1 -Az 2 



(6) 



In this edition of the tables the power of two is mistakenly written as j instead of 
n. In order to check the formula it is enough to expand the binomial into powers of z, 
integrate term by term and compare the coefficients with the formulae given in p5(| . 
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the function 



p(p,p) = T,Mp)Kr n (p') (7) 

n 

can be interpreted as a well-defined single-particle density matrix in the mo- 
mentum representation. The functions ip n are its eigenf unctions and the 
numbers \ n are the corresponding eigenvalues. Let us choose in particular 
for the functions ip n the energy eigenf unctions in the momentum represen- 
tation of a one-dimensional harmonic oscillator and for the numbers A n a 
geometric progression [El BJ: 



a 



v/?2 n n! 



cxp 



a 2 p 2 



H n (ap); X n 



:i-z)z r ' 



It is assumed here that \z\ < 1, so that the eigenvalues A n are correctly 
normalized. The parameter a is positive. Its interpretation in terms of the 
mechanical parameters of the oscillator is of no importance here. This choice 
corresponds to a single particle density matrix which using (0J) can be written 
in the form 



P(P,P') 



a(l — z) 



i« 2 (p 2 +p' 2 ) 



e 2 



$(ap,ap',~), 



(9) 



'71 T 

or substituting the definition of the function $ and replacing the parameters 
p, p' by the parameters K, q according to (0) 



p{K,q) = -7=\/t- -^exp 
% V 1 + z 



1 — z a 2 q 2 1 + z 



-a 2 K- 



l + z 



4 1 - z 



(10) 



This density matrix differs only by notation from the density matrix (0). The 
correspondence between the parameters is 



1 



a 



v + 1 



Following we have introduced the notation 



2AR > 1, 



(12) 



where the inequality follows from the uncertainty principle. 

4. In order to calculate the function G n (p, p 1 ) let us note first that it is the 
p, p' matrix element of the n-th power of the matrix p. Therefore, it has the 
same eigenf unctions and its eigenvalues are n-th powers of the corresponding 
eigenvalues of matrix p. Thus 
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k=0 



a(1 " z)W e-^ + ^<&(ap,aj/,^ 

'7T 2 



(13) 



Substituting the definition of function $ and changing variables from p, p' to 
i-C, q we find 



a(l — z] 



: exp 



-a 2 K 2 



1 - 

1 + Z r ' 



1 



^(1 - z 

It is convenient to rewrite this formula in the form [O 

i^ 2 1 



2 2 1 + Z 

-a q 

4 y 1 



2'- 



(14) 



G n (K, q) 



exp 



2A 2 



By comparison with the previous form A n = a n A 2 ; i? 2 = a n R 2 with 

1 1 + z n l(v + l) n +(v- l) n 



and 



vl-z n v (v + l) n - (v - l) r 



2 n ^ 



(u + l) 2n - (v - l) 2n 
5. Let us add some comments. 



(15) 



(16) 



(17) 



As mentioned, the transition from the one-dimensional case to <i-dimensions 
is very simple. The density matrix (2) gets replaced by a product 



P(P,P') = n 



1 



2ttA 2 



exp 



2A? 2 



Consequently, the integrals G m from (1) factorize and one obtains for 
each G m a product of expressions (15), in general with different values 
of A 2 and v for each factor. 

For the harmonic oscillator the energy eigenfunctions in the momentum 
representation and in the coordinate representation have the same form. 
Besides the replacement of p by x, the interpretation of the parameter 
a changes, but this is irrelevant for the present analysis. Thus, the 
results given in the present paper can be used, after some changes in 
notation, also for the study of Bose-Einstein correlations in ordinary 
space. 
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Since the energy eigenvalues of the harmonic oscillator are equidistant, 
the corresponding Boltzmann factors form a geometrical progression. 
Thus, formula (9) as well as the corresponding formula in the coordinate 
representation, can be used to calculate the canonical density matrix 



for the harmonic oscillator. According to Landau and Lifszyc |L6 
this calculation for the diagonal elements of the density matrix in the 
coordinate representation has been first performed by F. Bloch in year 
1932. 

The author thanks Andrzej Bialas for collaboration and stimulating dis- 
cussions. 
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